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1. Introduction 
1.1. The problem 

Topological defects are defined algebraically as operators commuting with the left and 
right copies of the chiral algebra Jy , and in particular are invariant under diffeomorphisms, 

[L n ,X] = [L n ,X} = 0. (1.1) 

We are concerned in this note with the 4-point crossing relation in the presence of defect 
operators: for trivial defects it reduces to the standard Belavin-Polyakov-Zamolodchikov 
duality relation for the correlators of local 2d fields. The presence of a defect line inserted 
between two local operators modifies their operator product expansion (OPE), resulting 
in creation of defect fields, and one is interested in the computation of the corresponding 
OPE coefficients. A special case of this extended BPZ relation has been exploited in to 
derive a general formula for the relative coefficients of the OPE of local fields of integer spin 
in the rational non-diagonal theories. Here we restrict to the diagonal theories where the 
computation of the duality transformation in the presence of defects is a straightforward 
consequence of two basic identities in CFT: the pentagon identity for the quantum 6j 
symbols (the fusing matrices) and the Moore-Seiberg , [|J torus identity, an equation for 
the 1-point modular matrix. 

As a side result of this computation, extended to the non-rational c > 25 Virasoro 
theory (Liouville CFT), one obtains an explicit general expression for the expectation value 
of the 't Hooft loop operator in Liouville theory, defined as the dual to that of the Wilson 
loop operator. The formula essentially reproduces the recently proposed ad hoc expression 
H, || and thus confirms the assumed duality of the two operators. 

The effect of some defect lines as creating "disorder" fields when attached to the 
local fields was pointed out already in 0. It was later thoroughly analysed in |7] and 
in particular the precise conditions on the type of the defect fusion algebra leading to 
Kramer s-Wanniers type duality in the rational case were described; for more on the defect 
fields from a TFT point of view see M, M. A topological defect interpretation of the 



Wilson loop operator in the 2d rational theories has been discussed e.g. in [[10]], The 
construction of the loop operators in H, H does not refer to defects.0 



1 The present work was in an advanced stage when it was announced [12] that a parallel work 
on the defect interpretation of the construction in 0, || is under way. 
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1.2. Preliminaries on the topological defects 

Let us summarise some of the consequences of the definition studied in 0, ||. 
In a rational CFT with a set {X 3 j} of representations the solutions of (|1 . 1| ) read 

*« = £ E V^r- p(j ' j;a ' a,) ' ( L2 ^ 

where P^^'^ a ' a ) are projectors in the representation spaces (Vj <8> Vj)' 3 , and the sum is 
restricted to pairs allowed by the nonzero values of the given modular invariant 

matrix Zjj, taken with their multiplicity (3 = l,...Zjj] \& is a unitary matrix of size 
As in the computation of the cylinder partition function, which leads to non- 
negative integer matrix representations (NiM-reps) rij a b of the Verlinde algebra [ 373 1 , JH 



n i n i = ^2-^ij S n s , T rij = nj* , j EX 



s 



Sji ^ , ^ 

Su 



iex,a=i,...z u 



one constructs partition functions on the torus Z Xn \ XlX2 _ Xn _ 1 (or, on the cylinder 
Zb\xiX2...x n a)i inserting arbitrary number of defect operators. They are sesquilinear (re- 
spectively linear) combinations of the characters of the chiral algebra representations with 
non-negative integer coefficients Vg.^ Xnl Xn G X (resp 7ij. Xl „, Xn . a ). The case of two 
defects on the torus leads to an equation for the multiplicities analogous to (|1.3|) 

T/- • T/- • — \ A \f. . *3 \f. . 33 t>. . \r. 1 — 7 ■ ■ T V- ■ — T/„ 

h ,33 

q C _ (1.4) 

~ y. y_ V 37 ju(l,l;a,a') ,T,(l,l; a ,a') * 

w ^j; 1 ~~ c ,c - x y 

l, l, a, a' 

The classification of the topological defects amounts in the classification of the NIM-reps 
(p. ,4j) . In the sl(2) related cases it confirms the results of Ocneanu, visualised by his gen- 
eralised ADE diagrams, with vertices associated with the set of defects |T^]. Another 
distinguished set of non-negative integers Vu- yx z = N yx z is provided by the identity con- 
tribution of left and right characters in the torus partition function with three inserted 
defects, Z z \ yx . This set serves as structure constants of an associative, in general non- 
commutative, algebra, the fusion algebra of defects (=Ocneanu graph algebra), as it allows 
to compute the fusion of two defects 

XyX x = ^2N yx Z X z . 
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The action of the defects on the boundary states 
X x \a) = S ^2n ax c \c) , 



c 

C?.i;a.0) (1-6) 



, . * 



j,a,/3 V^i^lj 



introduces another set of non-negative integers {n ax c }, interpreted as the multiplicities 
n i-x-a C = of the identity character in the cylinder partition functions in the presence 
of one defect, Z b - X)a = Y^j nj ]X - a b Xj( T ) = *E,j(nxnj)a b Xj( T )- Combining relations ( |L5D and 
( p..6|) implies that this set of matrices provides NIM-reps of the Ocneanu algebra 

n x h y = s ^N xy z h z . (1.7) 

Thus the computation of the partition functions on the torus and the cylinder with an 
arbitrary number of defects inserted is reduced to the knowledge of several basic structure 
constants, or, equivalently, the knowledge of the sets of unitary matrices {ip,^} in (|1.3|) , 
(ED, e.g., 

V u 

In the field interpretation, these multiplicities encode the possible holomorphic - antiholo- 
morphic content of defect fields; in general they correspond to non-local 2d fields. 

In the case of a "diagonal" theory, i.e., described by a modular matrix with Zjj = 8ji, 
the sets of defects and boundaries can be identified with the set X of representations of the 
chiral algebra and n, h, N coincide with the Verlinde multiplicities A/", while 

Vij = MiMj , v m * a;/. (1-9) 

etc. Accordingly the matrices diagonalising these multiplicities reduce to the modular 
matrix, #i J ' j);Q ' /3 = 8fi6 ul 6p ±S X j, and the eigenvalue in the r.h.s. of (|1.2|) is expressed by 
the ratio of modular matrix elements S x j/S±j. In the diagonal case the action of the defect 
operator X x (on (Vj <8> Vj) coincides with the action of its chiral analog acting on Vj (one 
to one with the Ishibashi states) 

*i = £l*£l'''*>&*l (1.10) 
3 lj k 
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In the WZW case this operator can be identified [fl l] with a generalised Casimir operator 
|T(J(see also [17] for the Virasoro minimal models) giving precise meaning of the Wilson 



loop operator [|T0|. The twisted partition function in the diagonal case can be interpreted 



alternatively in terms of related operators, associated with the two cycles of the torus 
X x (a)Xp(?) = Tr p (e 2 ^°- c / 24 )xi) = ^ Xp (f) 



This reproduces the monodromy operators in the derivation of the Verlinde formula [18 



In the Liouville case the sum in (|1.2[) is replaced by an integral over a continuous series 



of representations. The modular matrices were computed in [Tj| in the two basic cases - 
with the second representation x also belonging to the continuous principal series, or, with 
x belonging to the c > 25 infinite discrete degenerate series. The latter is parametrised 
by a pair (to, n) of positive integers, i.e., x labels a factor representation of a degenerate 
Verma module of scaling dimension A(x) = x(Q — x), with x = x mn = — ^^-b — , 
Q = b + ^. Thus, there are two types of defects in the Liouville theory, with two different 
ratios of modular matrix elements in ( |1.2| ), corresponding to the two types FZZ |2(J or ZZ 
T^d , respectively, of boundary states; this has been discussed in detail recently in [^T| . 
For the purpose of the comparison with recent work [|5| , on the relation of Liouville 



theory to the 4d supersymmetric gauge theories ||22|| , we perform in more detail also the 
explicit computation in the Liouville case. We will restrict mostly to the ZZ case, a quasi- 
rational theory, similar in many respects to the rational c < 1 theory. The main section 
2 deals with the diagonal rational case. The technical difference with the Liouville theory 
is in the use of different normalisation of the chiral vertex operators - traditionally in the 
minimal models the Z2 symmetry is effectively fixed, the reflected operators are identified, 
and in particular one can use bases of conformal blocks for which the fusing matrices are 
unitary. The three parts of section 3 deal with the Liouville case and contain: a collection 
of basic formulae (sect. 3.1); more details on the Liouville torus identity and the crossing 
relation in the Liouville case (sect. 3.2); an explicit example compared with the proposed 
expression for the expectation value of the 't Hooft operator in 0, (sect. 3.3). The 
Appendix recalls the OPE formula of P] . 



2. The duality relation in the rational case 



For simplicity of notation we shall restrict to the sl(2) case, but will keep the conju- 
gation so that the higher rank generalisation is straightforward. In the rational case one 
can choose a unitary fusing matrix F. 

As in || we consider a 4-point function of four local fields with insertion of two defects; 
for the purpose here we shall restrict to the diagonal case with scalar fields; the second 
label will be suppressed, 

G = (0|$ a4 (5 4 ,l 4 ) $ (a3 (5 3 ,l 3 ) X x $ a2 (z 2 ,h) $ ai (zi,2i) X x \0) 

^^l^( a 4,a 3 ,a 2 ,ai;2)| 
, S ij S n (2.1) 

£ 7 (a 3 , a 2 , ai, a 4 ; z)Qs(a s , a 2 , a\, a 4 ; z) . 



■iryy s : j i- 



*-i ^ " 

7,5 3 



a 4 ai 

a 3 a 2 



771* 

*j6 



a 4 ai 

°3 a 2 



In the second line we have used that in this channel the defects are diagonalised, applying 
the OPE expansion of the local fields and inserting the diagonal version of the general 
formula for the 2-point function 



(0|$ (J * ia) X x $ (J , i/3) |0) = 5 j;j ,8 



(J;a,/3) 



3,3 \f,J 



l$(J-,a)$(J,/3)|0> 



(2.2) 



In the third line of ( |2.1| ) we have performed the braiding of the chiral blocks. In the new 
channel the defect modifies the OPE expansions of two local fields 



G = (0|$( as (z 3j z 3 ) X x $ a2 (z 2 , z 2 ) $ ai (zx, zt) X x $ a4 (z 4 , z 4 )\0) 
It remains to perform the initial summation over j in (|2.1|), i.e., 



A 7,S 



^ Su tjl 

3 J 



a 4 ai 



a 4 ai 
a 3 a 2 J 



, with A^ } s = 5 lS . 



(2.3) 



(2.4) 



The first step is, using a standard relation derived from the pentagon identity, to rewrite 



the product of F matrices (|2.4|) as 



F 

r 31 


CL4 CLl 


T7T* 

*jS 


a 4 ai 






_ a 3 a 2. 




_ a 3 a 2. 





da\ d a2 d a3 d a4 



aiaZ 



a 2 7 
j at 



a 4 a,2 



«3 5* 
j ai 



(2.5) 
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We have also used the relation derived from the pentagon identity, taking into account the 



symmetries of the 6j symbols and using that Fn 

F\rn 



a a 
a a 



l/d a : 



1 i ' 




~i* f 


F m l 








.3 3. 



dm J\J"..m 

didj Jl 



(2.6) 



It implies that for a unitary F the values of Fi c and F y i are determined up to a sign in 
terms of a square root of a ratio of q-dimensions; the positive sign is chosen. We next 



apply the pentagon identity itself, representing the product in (|2.5| ) as a sum 



«3 S* 

3 ai 



7* 5 ' 
a 2 a 2 



a 2 7 
3 at 



F * 



7 5* 
ai ai 



y j 



a 2 7 
a 2 S* 

a 2 ai 
a 2 al 



a 4 y 



F a *i 



" 7 (5*' 


F 

± a 4 a 2 


a 2 


y ' 


ai ai 


. 3 


a 1 _ 



3 3 

CL2 CL2 



(2.7) 



a 2 a>2 
a 2 a 2 



Inserting this into Q2.4Q we can now perform the summation over j, using the MS torus 
identity M, 0]. It can be written e.g. as 



s 

= S qi (p)F rp 



qs 



31 J2 

32 31 
q q 



e 2 «( A i- A i.+ A ii+Aj 2 -A m )^ 



* 32 

iji 



mp 



3i 32 

i i 



[2.1 



Taking r = 1, hence s = 1, gives an expression for the modular matrix S(p) of one point 
functions on the torus 



•lp 



3 3 
3 3 



e »T(2(A,4A r A m )-A p )^ if 



3 * 
3 i* 



3 3 



(2.9) 



and Sj^p) = e l7T p Sj*i(p*). Inserting the transposed version of ( |2.9| ) in ( |2"T8| ) - taken for 
p = 1, and inverting with F^, 1 gives a "fusion" like representation for the product of two 
S-matrices, 



Fy*l 



dx \ S x j p 

a.2 a*] ^-f Su 

a 2 a,2 



y 3 



a 2 ai 
a 2 a\ 



F a *i 



3 3 
a 2 a 2 



Sa 2 x(y*)_ e i 7T /\. u S aiX (y) 



Su 



(2.10) 



In particular, for y = 1 (|2.10[ ) reproduces, using (|2.6|), the Verlinde formula for the fusion 
multiplicity Af. Combining ( |2.5|) , ( p.7|) , (|2.10|) we get finally for 



d x A 



(x) 
7,5 



d-ydfi 



da\ da,2 da,3 d a4 



/ j a 3 y 



7 6* 



S*a*Ms aiX (y) 



Su 



Su 



7 5* 



(2.11) 
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The range of y in (|2.11|) is determined by the multiplicities M xy x ,M ay a of the 1-point 
blocks on the torus, transformed by the 1-point modular matrix S ax (y), a = ax, a^. We 
see that we can interpret y as a defect label, appearing in the defect fusion multiplicity 
M xy x . In turn the values of y restrict the possible pairs (7,5) of representations labelling 
the left and right conformal blocks in ( |2.1|) as dictated by the multiplicity A/" 7 <5* y , implicit 
in the F matrices. Altogether this corresponds to the first of the relations in fll.8|) for 
z = x, 



y y 



(2.12) 



In other words, Vj 1 s*;i v is the multiplicity of a defect field vQ^^z, z) that is created in 
the OPE of two local operators, modified by the inserted defect line operator X x as in 



— A" 1 _-A s 



-,, E d »i^» (7 '' )iWW *(7,«)(^^) + "- (2.13) 



and the OPE coefficients in (|2.13| ) are determined from ([2.1 1|) up to a normalisation of the 
defect field 2-point function. Alternatively, inserting 



~x5*~ 


t y'i 


~x5*~ 






x 7 




x 7 



we can rewrite (|2.11| ) in terms of the modular matrix transforming a 2-point chiral block 
on the torus 



Sw, g (7i,72) = F~\S® 1)F = J2 F y<i 



Xl2 

l x 7i J 



Si X {y)F ai 



7i 72 
i i 



(2.14) 



In the higher rank cases it depends on four more indices, their range being determined 
by the fusion multiplicities M lia l ,N l2 i a and M liq x ,N l2X q of the chiral vertex operators 
involved. For ( |2.11|) we have 



d x AW = y/djr s J2 



S a*,a 3 ;x,qhi S*) S ai >a * ^(7, 5*) 



,, \/da 2 da 3 \fd ai d a4 

The summation here corresponds to (|2.12|) rewritten as 



1 

5*x 



(2.15) 



(2.16) 



Taking 7 = 1 = S in ( |2.11|) represents the leading contribution of the identity block in 
( |2.1|) , with the restrictions 03 = a 2 , 04 = a*. Using ( |2.2| ) this reproduces indeed the r.h.s. 
of ( |2.11| ) for these values. 



In the analog of the chiral interpretation ( |1 . 1 1| ) here, the modular transformation of 
the characters is replaced by braiding (fusing) of the conformal blocks on the sphere 



X x (a)G p (z) = (0|^ >0 .(zi)^Sp(^)-X'X > a 1 (^Xo(2i)|0> = G p (z) 



(2.17) 



s 3 



and formula ( |2.11|) gives alternative expression for A 



(x) 

ps ■ 



The MS torus identity has been encountered in the boundary CFT in |23], [14]; it 
has been observed that in the diagonal case the two Cardy-Lewellen bulk-boundary eqs 
|24j], J25| both originate in this identity. The first equation corresponds to ( |2.8|) (with 



boundary labels q, r and F rp , F qs substituted by the 3j symbols - the boundary field OPE 
coefficients), while the second equation is identified with the (transposed version of the) 
relation ( |2.10| ). Accordingly the bulk-boundary structure constant in the diagonal theory 
is proportional to the 1-point modular matrix. 

The duality transformation relating the correlators (|2.1[ ) and ( |2.3| ) has been recently 



discussed in following and comparing with the permutation brane approach; this 



consideration does not yield, however, explicit formulae like (|2.11|), (|2.15p. 



3. The Liouville case 

3.1. Collection of Liouville formulae 



The quantum 6j symbols F for the Liouville theory have been computed in [£6] , [27 
It is convenient to change the normalisation of the chiral blocks, ^ F 'Gp — > Gp 



(F) £/3(«4, o> 3 , «2, cti]z) = N(al, a 3 , 0)N(0, a 2 , «i) Gp{ot±, a 3 , a 2 , a x ; 

T b (Q)T b (2/3 1 )T b (2^ 2 )T(2Q - 2#$) 



N(0 3 ,0 2 ,0! 



r b (2Q - /3 12 3)r b (/3 1 3 2 )r 6 (/3 2 1 3 )r 6 (/? 1 2 3 ) 



(0123 = 0% 1 0f 2 = + — 03, etc.) so that Qp transform with the matrices 



G 



Ps ,Pe 



03 02 

04 01 



N(^0 3 ,0 2 )N(0 4 ,0 6 ,0 1 ) 
N(0 4 ,/3 3 ,0 5 )N(c H ,0 5 ,0 1 ) 



Ps ,p6 



03 02 

04 01 



(3.1) 



(3.2) 
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The Liouville bulk 3-point constant is then given by 

-Q 



N(/3 3 , fa, fa)N(Q -fa,Q-fa,Q-fa) = 2n\— J] W(Q - fa) C{fa, fa, fa)' 1 • (3.3) 



Here A := tt^ T(b 2 )/T(l - b 2 ) b 2 ~ 2b and we have used the ZZ variable p 



^> = i^^A^(=-W(iP)«). (3.4) 



Recall that the product of W(a) in ( p.4j ) and its reflected counterpart W(Q — ct) is pro- 
portional to a modular matrix element, while the ratio gives the bulk reflection amplitude, 

W(a)W(Q-a) = J,^^ = -4sin7r6(2a-Q) sinf (2a- Q) =: S 0a 

W(Q-a) _ T b (2a) _ l ' J 

W» T,(2a - Q) W ' 

More generally, for the case of a degenerate representation x and generic charge a the 
modular matrix reads |0| (up to an overall normalisation) 

S Xm , n a = -4sin7r6m(2a-Q)sin^(2a-Q) = S Xm nCi - S x _ m ^ a (3.6) 

where 

5^ = 2cos7r(2a - Q)(2/3 - Q) (3.7) 

is the FZZ type modular matrix, computed for two generic representations. 

The Weyl reflected charge in the second line of ( |3.6|) corresponds to the only singular 
vector at generic b 2 of the reducible Virasoro module of highest weight A(x m>n ). This 
relation, coming from the character formula for the degenerate representations, extends 
to other quantities of the theory, e.g., the corresponding fusion multiplicities, Map 1 and 
M 7 

M ax 7 =A r QX J -Af a x 7 • (3.8) 

Here the l.h.s. is a finite sum of delta functions, while Map 1 is given |TI| by an integral 
formula of Verlinde type, i.e., it is diagonalised by y/2S a s in ( |3.7|) and its eigenvalues 
(1-dimensional representations) are given by the ratios S a s/Sos. 

The F matrix is invariant under reflection fa — > Q — fa of any of the indices ||26|1 , 



equivalent to a complex conjugation for pure imaginary iP = Q—2fa Extended to arbitrary 
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values of the charges, the gauged G matrix ( |3.2|) satisfies the standard symmetry relations 
with the star operation understood as a reflection, j3* = Q — f3 



G 



Ps A 



p3 

Pa 



02 
Pi 



G 



p 5 ,pi 



Pt 
Pt 



Pi 

P2 



G 



PS ,06 



P2 p 3 
PI PI 



(3.9) 



The locality of the scalar 4-point function is rewritten in terms of the G matrices as 



S. 



70 



«4 a 1 

CK3 Q.2 



G 



Pi' 



CK4 OL\ 
CK3 OL 2 



5(7-7% 



(3.10) 



The integrals here and below run along Q- + ilBL + . We shall exploit the relation of the fusion 
matrices to the Liouville boundary field OPE coefficients C [BSL 



G 



T2 ,Q~Ps 



P2 
03 



Pi 
01 



N{Q - P 3 ,P 2 >Pl)R(<T3,Q ~ P3,<Tl) 



R(a 3 ,p 2 ,a 2 )R(a 2 ,p 1 ,a 1 ) 
where R is the ratio of the two gauge factors; 



C, 



<?2 ,Q — @3 



P2 



Pi 



R 1 (a 2 ,7, cr 3 ) 



£(0-2,7,03) 



A 



■>+"3-"2 S b (~f + a 2 - 03) (7 + 03 - 02) 



(3.11) 



(3.12) 



^(02,7,03) ' " 5 b (2 7 ) 
The OPE coefficients are related to the coefficients of the boundary field 3-point functions 



/"fO"3,°"2,0"l 
U P 3 ,P 2 ,Pl 



5(0-3, P, 0l) 



a 



<?2 ,Q-P 3 



P2 Pi 

03 01 
0(03, Q - P, 01) 



5(0-3, P3, 0l)C<T 2 ,p 3 



P2 
03 



/3i 
01 



(3.13) 



g(?3,P, 01) 

with the boundary reflection amplitude p0] defined in the second line. In the case when 
the three charges Pi in ( |3.11|) are constrained by a charge conservation condition, the 3- 
point function Cp^^'pi develops poles. The residue corresponds to the correlator, which 
can be computed in the half-plane Coulomb gas formulation of |2C]. We shall denote it 
and the residues of the corresponding G in ( |3.11| ) by the same letters. For ^7 Pi — Q = 
(absence of screening charges) the residue is 1, so in these cases G reduces to the gauge 
factor in ( |3.11|) . Furthermore any C related by a reflection to a trivial one is also simple, 
being obtained by applying the boundary reflection matrix as in (|3.13| ). This modifies one 
(or two, or three) of the ratios ( |3.12j ) in ( |3.11| ), replacing 0(0-4, 7, 03) with 0(04, Q — 7, 03). 
Examples of G matrix elements obtained this way will be used below: 



G 



PQ 



G 



07* 



or a 

7 7. 
'a* (3 

_Q 7 



1 



simibQ sin j-Q 



G 



«7 



sin7ro(2a - Q) sin f (2a - Q) 
7" 



G 



po 



7 7 
a a 



a (3 



(3.14) 



1 
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where the quantum dimension d a = appears; furthermore 

sin7r6 2 



G 



Q+b,a±b/2 



Q 



G 



a±b/2,Q+b 

More generally, denoting 

S{a 3 ,p, at) 



Q+l 



a 



_b 
2 
a 



sin7r6(2a — Q) 
sin7r6(2a =p Q - Q) 



(3.15) 



sm7r 



2b 2 



G 2 (cr 3 , P, a t ) 



S b (2(3 - Q) 



S b ((J 2 + &1 ~ (3)S b (Q - /3 + (T 2 - (T\) 



(3.16) 



W{Q -P) ^' S b (P + a 2 + a 1 - Q)S b (P + a 2 - a x ] 

we can write a compact formula for the general Coulomb gas boundary coefficients C, 
obtained as a residue from the Ponsot-Teschner (PT) formula P^j , 

Pi Pi 

S b (2p 2 + mb+ m(2ft) ™ " G 2 (a 3 - 



a 



°2 ,Q~f3 3 



27r Res / 3 1 



r (PT) 
-Q+mb+n/b=0 ,Q-/3 3 



Pi Pi 
03 C\ 



g(Q-P 3 ,P2,Pi 



S b (2p 2 )S b {2p x +mb + 



EE 

b > k=0p=0 



(k—m)b 



2b 



;Q-Ps^ 



(k — m)b , p—n 



2b 



G 2 (a 3 -f-^ 5 ,Q-p 2 -^-^ E ,a 2 ) 



Gi(cr 3 , Q - Ps , cri] 

' _l\m(p+l)+n(k+l)+r 



G 2 (CT3:Q- Pl,Vl) 



S b ((k+l)b)S b ((m-k+l)b)S b (mS b 



■ b 



(3.17) 

Further some of the Pi charges in (|3.17|) can be set to degenerate values; it is a polynomial 
in the boundary parameters 2cos7r&(2o"i — Q) ,2 cos ^(20^ — Q). Rewritten in terms of 
finite products of sine-functions ( |3.17[ ) admits analytic continuation to the region c < 1 
|29H and in this sense the integral formulae of [p6|, |28| are universal. 



The following relations follow from the pentagon identity, 



G r 



' j h~ 


— G b k* 


~i* j' 


b a 




a c 




= Gi c * 


~k* a" 




.J b 



GcQ 


~ k* k 
a a 


GbQ 


_ a aj 


di 
d c 





G 



b*k* 



J i 
c* a* 



di_ 
dk 



(3.18) 



where in the second equality we have used the (Coulomb gas) values ( p.!4[ ), particular for 
the chosen gauge, and the third equality is obtained repeating the first one. This relation 
is derived alternatively by using ( |3.11|) and the cyclic symmetry of the boundary 3-point 
coefficients in the l.h.s. of ( |3.13| ). In particular ( |3.18| ) implies 



Goi 



j k 
jk* 



G 



OQ 



k* k 
k* k* 



G jQ Vk* k*] 



G 



Qi 



k j 
kj* 



dj_ 

dk 



(3.19) 
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Here we have replaced delta function singularities on both sides with the residue values; 



for the precise details of treatment of these singularities see Appendix B of [[28]. 

In all these relations it is assumed that the triples of representations are consistent 
with the corresponding fusion multiplicities. With this data the analogs of the first two 
steps (P2D> (P- u) in the rational case now read 



G 



Pi 



Ct4 Cti 
«3 Oi 2 



G 



06 



d^ds 



dyG a * y 



(X4 Cti 

a$ a 2 
"7* 5 
a 2 a 2 



Gctia^ 
Galy 



a 2 7 
(3 ol\ 
7 



G 



y*0 



a% 5* 
(3 ai 

«2 Oil 

a 2 a\ A 



■j as 



dots ^C*2 



G 



a*Q 



(5 0* 
a 2 a 2 



(3.20) 



3.2. The torus identity and its application 

The basic MS torus identity in the Liouville theory is an integral relation 

S rx (s) J dme 2 ^ A{x) - A{m)) G & 

__ f> in(A(p)-A(j 1 )-A(j 2 )) 



XJ2 



G 



rap 



3\32 
X X 



dqS qx (p)G 



sq 



rj2 
rji 



G 



rp 



32 3i 
q q 



(3.21) 



The identity is gauge invariant and can be rewritten in terms of the F matrix with 

{F) S ax (s) = S cex (s) Z[ a,S,a ] = (F) S a AQ - s) • (3.22) 
N{x, s, x) 

It depends on the range of the representations, here symbolically written in general form; in 
the FZZ case the notation S a p(p) will be used. Setting r = = s, hence q = 32= j\ = «, 
one obtains expressions for the 1-point modular matrices, integral in the generic case. In 
the case of our main interest, when x (and hence p) in (|3.21|) is degenerate x = x m ,n, the 



first integral in (|3.21|) is replaced by a finite sum with the Coulomb gas expressions of the 
G matrices appearing, 



Sax(p) 



Sqx 



G 0p [ QQ *1 ^ 
Sjo 



^ ^(2(A(x>hA( a )-A(u))-A(p)) G ^ 



a x ' 


n 

^up* 


' X* 


X 








.a x* . 


a 


a 



o j0 \^ i 7 r(2(A(x)+A(a)-A(u))-A(p)) r f 

P^i IX* X* J u 



a x ' 


G u o 








.a x* - 




a a 



(3.23) 



(where u = a + kb + l/b with the range of the integers fc, I restricted by the fusion rule 
of the degenerate x m , n )- The second equality follows from the pentagon identities. Note 



that S ax (Q) = ^-S ax from ( |3.22|) . As in the rational case, (|3.23|) provides for p = an 
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alternative formula for the modular matrix ( |3.6|) , which can be also written in integral 
form 

So Im ,„ = SOO / rf7e 2-(A(a) + A(x m ,„)-A( 7) ) ^^^7 (3.34) 



This and the analogous formula for S a p/ 'Sqq with M olXm/n 1 replaced by A/" a /3 7 can be checked 
as in the rational case. In the generic integral analog of ( |3.23| ) S ce /3(p) stands in the l.h.s. 
of ( 3.23 ) and an alternative representation can be obtained similarly to the computation 
in [|30| of the FZZ bulk-boundary constant, by solving the set of finite difference equations 
obtained when setting j 2 = —6/2 in the integral analog of ( |3.23| ), equivalent to the basic 



identity ( |3.21|) ; see also below. 

The expression ( |3.23|) simplifies if 2a = p*, or p (or if 2x = p, or p*) since, as 
discussed above, the G matrix in the r.h.s is simple, does not involve a summation and 
reduces to a product of gauge factors. The simplest examples are provided taking ( |3.23| ) 
for x = £2,1 = —b/2 and p = —6, or p = Q + b. For these values the sum contains two 
terms, u = a ±6/2, e - i7rA (- fe ) = _ e nr2Qb anc j one CO mputes the 1-point modular matrices 
using the fundamental G matrices ( |3.15| ), ( |3.14| ), ( |3.19| ) 



S h ' _ d a 2ie inQb sin7r26a sin7r6(2a - 2Q) 



Sqo 

S b (Q + b) 

01 ~2 



d. 



sin 7T 



b' 2 



S 



00 



d_b 
2 



e inQb 2i simrb 2 . 



Next we set p = in ( |3.21| ) and apply ( |3.23|) for the sum in the l.h.s. of ( [3.21 
also ( |3.14| ), or, changing notation, r = a.\ , a = a 2 , q = f3, we get 



d/3 -^-G y *p 



S, 



00 



Ct2 Oil 

ai2 a* 



G 



a.2 OL2 



S a2X (y*) e l7rA(y) S aix (y) 



S, 



00 



dl 



S 



(3.25) 



using 



(3.26) 



00 



This relation is what we need when evaluating the analog of (|2.4j ). Indeed we have an 

additional factor Spo, coming from the measure in ( |3.10| ), which cancels the denominator 

s 

of the defect eigenvalue Combining with ( |3.20|) we finally obtain 



A 



(x) 
7,<5 



d/3 ^ Gp-y 



S 



00 



«4 «i 

a% «2 



G 



06 



«4 a± 

Q!2 



djds ^ ] G 



a,y' 



7* 5 
a.2 ol 2 



S a2X (y*) e^ A ^ S aix (y) 



Soa 



d. 



Soa 



G 



7 5* 



(3.27) 
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We have replaced the integral in ( |3.20|) by a sum once again using the fact that for degen- 
erate representations the fusion matrices are represented by residues of the initial singular 
expressions. Using that S ax (Q) = j^Sax and (|3.14|) one checks that for trivial defect x = 



(hence y = 0) (|3.27|) reduces to = d 7 <5(7 — 5) in agreement with ( |3.10| ). 



For x in the continuous series, the integral analog of (|3.27|) holds, with the ratio 
in the l.h.s. (confirmed as a defect eigenvalue), while S a2X (y*)/ Sq x ,S aiX (y)/So x will 
appear in the r.h.s.; the remaining q-dimension factors are unchanged. 

The generic analog of ( |3.26| ) considered as an expression for S a2X (y*) simplifies for 
the choice a.\ = y/2 of the other charge by the mechanism discussed above: this is more 
transparent in the transposed version of ( |3.26| ) obtained using the identities ( |3.18|) and 
S ax (v) = tSzctv*). On the r.h.s. appears the constant Res 2ai=y S xctl (y*), while with 
this choice the G matrix in the l.h.s. of (|3.26[) is replaced by a Coulomb gas correlator 
from (pip : 



G 



OL 2 



a 2 Q—% 



Q 2 



P (3* 
a 2 a 2 



(In 



Q 



y y 

2 2 

a 2 a 2 



Soo S h {f3 + Z-a 2 )S b {(3+l+a 2 -Q) 
S 2 b (y) S b (f3 - | + a 2 )S b (P - | - a 2 + Q) 



(3.28) 



The transposed version of ( |3.26|) (taken with either of the two modular matrices S x p or 
S x p) is to be compared, with the corresponding Cardy-Lewellen type equation for the bulk- 



boundary reflection coefficient. Special cases of this equation have been used in [P0| , |19 



for the determination of the half-plane bulk 1-point functions in the two Liouville cases. 
The general equation has been exploited in |31| to give an alternative derivation of the 
FZZ bulk-boundary constant R x (a,y), computed in |30f| . In obtaining (|3.28|) here we have 
followed a similar argument to that given in |31]] - indeed one recognises in the expression 
( |3.28| ) the Fourier transform R(a 2 , y; ft — Q/2) up to /3-independent factors. The resulting 
expression for S xa (y) can be written in terms of 6-deformed hypergeometric functions (see 

B H) 



S xa (y)e^ A{y)/2 = ^S ax (y*)e^yy 2 



5 ^^^QF ei7r(±(2Q ^^ )(2a ^ )jF6(y ' ^-^^ M2a-Q)+Q; 2x-Q) 
S %t-Q) ) ( ^ e^+y-^~^F b (y, 2a + y-Q; 2a; ±(2x - Q)) . 



(3.29) 
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For y = it reproduces S xa in and the residue of S xct (y*) at a = y/2 is consistent 
with the FZZ analog of eqn ( |3.26|) at this value. The required symmetries like the one in 



the first line of (|3.29|) are checked exploiting relations for the 6-deformed hypergeometric 
functions. Note that the ratio 



s xa (y)e l7rA(y)/2 



tA( V )/2 



W(a)g(a,y,ai) W(a)g(x,y,x) 
has the correct properties under reflections with respect to the bulk a and the boundary 
y field charges as required for the FZZ bulk-boundary constant R x (a,y). 

For (odd) degenerate y = y2m+i,2n+i = —mb — n/b, one furthermore obtains, account- 
ing for the residues of the poles of the integrand defining the hypergeometric functions a 
formula for degenerate x = x rjS : the computation is analogous to the derivation of the ZZ 
bulk-boundary coefficient in [[31]. Namely, denoting 7^/ = kb + l/b, 

dry 



( (y) = S ax (y*) 



S Xr a{y) S x 



[-l) kn+ml S b ((m + l)6)5 b (2±i; 



S b (2a) 



m n 

_ i-rr A(y) v — ^ ^ ^ 

f: 2 > > S 

k=Q 1=0 



V. X 

2' 



S b (2Q-2a) 



S b ((m + 1 - k)b)S b {Z±£=l)S b {(k + l)b)S b (^) S b (2a+ lk>l ) S b (2Q-y- lKl -2a) ' 

(3.30) 

The values of y are restricted by the degenerate fusion multiplicity, < m < r— 1, < n < s— 1. 
The formula ( |3.30| ) provides alternative to ( |3.23| ) representation for S aXr s (Q— y2m+i,2n+i)- 



It is checked to reproduce the particular example in (|3. 25 



3.3. The 't Hooft operator: example 

Here we compute ( |3.27| ) for the simplest example considered in ||, 0, namely x = 
—6/2, so that y takes the values y = 0,-6. 

The G matrices for 5 = 7 ± b are all straightforward to compute, being related as 
explained above, to trivial boundary OPE coefficients, 



G 



03 03 



G 



d_ b /2<l a3 sin7r2a"3& sin7r6(2o"3 — 2Q) 
d^+t, sin 71276 sin 716(27 — Q) 
d_ b / 2 d- b sin 716(7 — a 3 + a 4 ) sin 716(7 — 04 + 03) 



03 0"3 



d 



-y+b 



sin 71276 sin 716(27 — Q) 



G 



a 4 ,Q+b 



Q—i 7—6 
03 03 



d_ b / 2 d a3 sin7T2a3& sin7r6(2a3 — 2Q) 
d 7 +b sin 716(27 — 2Q) sin7r6(27 — Q) a4 

sin 716(03 + 04 — 7) sin7r6(2Q — 7 — 



Q-7 7-& 

03 0"3 



d b d- b 



04 - 0"3; 



sin7r6(27-2Q) sin7r6(27-Q) 



(3.31) 
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For the only matrix that involves two terms (the case m = 1 , n = in ( j3.17|) ) one gets 

d-b cos 71-6(27 — Q) cos 716(203 — Q) + cos Tib{2o± — Q) cos7r6 2 



G 



0-4,-6 



Q-l 7 

03 0"3 



sin 7t2<7 3 6 sin 716(203 — 2Q) 



d 



7 7 



(3.32) 



Besides ( |3.25| ) we also need 
S b 

2 cos7r6(2c>!i — Q) 



S b (Q) d b S 



S, 



00 



da 2 



d 6 2cos7r6(2a 2 -Q) . (3.33) 



Let us change the notation (a 3 , a 2 , «i, 0^4) — >■ (04, 03, (T 2 , 01), so that B^\{y), as defined 
in ( [3.271) , reads 



7* 5 
C3 03 



5. 



00 



S, 



00 



7 5* 

C 2 CT 2 



(3.34) 



Collecting all formulae obtained from ( |3.31|) , ( |3.32| ) with the proper conjugations and 
change of variables, and using also (|3.14| ), ( |3.25| ), ( |3.33| ), we obtain for ( |3.34| ): 



d(-6/2)/ ,x _ 4sin7r6((2 + o-3 i -7) sin7r6((5+cr|-7) sin7r6(2(5-7-cri 2 ) sin7r6((7i 2 -7) 
7 > 7 " b ^~ >~ sin 716(27 - Q) sin7r6(27 - 2Q) 



,(-6/2)/ ^ _ 4sin7r6(Q+7-a34) sin7r6(7+a 3 4-(5) sin 716(7+0-2) sin7r6(7+o 2 ) 



sin7r6(27 — Q) sin7r627 



^(-b/2)/ 7\ COS7r6(27 — Q) COS 71-6(2(73 — Q) — COS7T&(2(J4 — Q) COS7T6 2 

7,7 ^~ ' ~ d_ k sin 7r62 7 sin 716(27- 2Q) 

2 

4(cos 716(27— Q) cos7r6(2o 2 — Q) — cos7r6(2oi — Q) cos7t6^ 



R (-b/2) 4 cos nb ( 2a 3 - Q) cos ^H 2ar 2 - Q) , 
fl w (°) = 1~, ' d 



(3.35) 



(3.36) 



; -6 = — 2 COS 7f6 . 

«_ b ~2 
2 

where 012 = oi+o 2 , o 2 = 01 — er 2 , etc. The above, normalised by <i_ b / 2 , expressions for 
By t s(—p)/d_b/2 should be compared with formulae (5.32-34) of ||. Apart from, presum- 
ably, a sign typo in (5.34): 

sin7r6(a — mi 2 — 6) shi7T6(a — 77734 — 6) — >■ sin7r6(a — mi 2 +6) shi7r6(a — 77134+6) 
the formulae coincide for (7; 04, 03, cr 2 , 01) — > (a; 777.4, "7,3, "72, Q— mi), i.e., up to a reflection 
of one of the charges. Here o\ = a 4 is the charge of the first vertex operator in the Wilson 



loop channel, cf. (|3.1| ). 
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The sum of the two terms in (|3.36| ) can be cast into a form which makes explicit the 
symmetry (02,03) — ► (o" 4 ,o"i). It coincides (up to a relative sign) with the term H in 
(5.39) of [j5| upon identification of the charges: 2a j — Q = 2irrij , 27 — Q = 2iP. In the 
initial basis of conformal blocks the duality relation reads 



d/3C(a 4 , a 3 , /3)C(/3*, a 2 , "Ott^ l (F) £/3(a4, «3, «2, u\\ z)\ 2 



I 

= fdyj d5j2C(a s ,a 2 ,( 1 ,6))C(( 1 \8n,a 1 ,a A )^^B^}(y)x (3-37) 
(F) £/ 7 (a 3 , a 2 , oti, a 4 ; z) (F) £|(a 3 , a 2 , «i, «4; 2) 



where we have denoted (consistent with ( |3.3| ) for 7 = 5) 



C(a 3 , a 2 , (7, <5))C((7*, 5*), ai, a 4 ) 
^2x^2 rr 4 



W( 7 ) 



(2tt) 2 A^ ni=i ^(<5 " «i) S. 



70 



iV(o!3, 7, a2)N(y, a 4 , ai)A^*(ag, 5, a 2 )iV*(5, a 4 , «i) 

iV(o!3, Q - 7, Q - a 2 )N(ai, Q - a 4 , 7) 



C(a 3 , a 2 , l)C{Q - 7, «i, ^4) 



iV(a 3 , Q — 8,Q - a 2 )N(a 1 , Q - a 4 , 5) 



(3.38) 



Multiplying ^1 7 ±6 ( — ^) in ( |3.35| ) with the ratio of N- factors, relative to the diagonal 
constant in the last line in (|3.38| ) , we get expressions invariant under any reflection cuj — > 
Q — a.{ of the four charges. These normalised expressions coincide with H± in (5.38) of 
(up to an overall factor 2ix) under the above identification of the charges. 



Appendix A. The defects and the OPE coefficients of local fields 

In the non-diagonal rational cases the identity contribution in the duality relation for 
the correlators (|2.1| ),( |2~T3"| ) is nontrivial and implies an explicit formula H for the relative 
OPE coefficients of local fields $j ;a (^, z) , I = (i, i) , a = 1, of arbitrary integer spin 

Namely (restricting to the sl(2) case) one obtains 



V d [K ;7) d (K;7 } — - — — (\ 9) 

fc,fc,7,7' 
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Using the unitarity of $ one gets an expression for the product of OPE coefficients which 
involves a summation over the complete set of defects for the given modular invariant .H 



The ^-ratios in (|A.2|) serve as 1-dimensional representations of an associative, commuta- 
tive algebra, dual to the fusion algebra of defects. This universal algebra generalises the 
Pasquier algebra ]33|] associated with each of the ADE NiM-reps in ( |1.3|) , which determines 
the subset of OPE coefficients with scalar labels only P3[ . 
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2 We use the opportunity to correct some inaccuracies in [01: formula ( |A.2| ) is slightly more 
general than what stated in [0J and reproduces also some of the signs in the case with non- 
commutative N, the Deven series: note that there are different bases for the pair of doubled fields 
in this case and accordingly different bases for the \& matrices. 
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